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1. Introduction 

Astrophysicists distinguish between three different types of compact stellar objects. These are 
white dwarfs, neutron stars, and black holes. The latter constitute a region of space which has so 
much mass-energy concentrated in it that no particles (not even light) inside the black hole's event 
horizon can escape the black hole's gravitational pull. The situation is very different for neutron 
stars and white dwarfs, which are about as massive as the sun (mass Mq = 2 x 10^" kg) but whose 
radii are much smaller than the sun's radius {Rq = 7 x 10^ km). Model calculations predict that 
the matter in the cores of neutron stars is compressed to densities ranging from a few times the 
density of an atomic nucleus, 2.5 x 10'*^ g/cm^, to densities that may be ten to twenty times higher 
^, ^]. In comparison to that, white dwarfs of average mass, M ~ O.6M0 are at least by a factor 
of 10^ less dense than neutron stars. The tremendous densities (and thus pressures) existing inside 
neutron stars make them superb astrophysical laboratories for a wide range of fascinating physical 
studies [[1], ^, ^]. These include the exploration of nuclear processes in an environment extremely 
rich of electrons and neutrons, and the formation of new states of matter, like quark matter which is 
being sought at the most powerful terrestrial particle colliders. If quark matter exists in the cores of 
neutron stars, it will be a color superconductor whose condensation pattern has been shown to be 
very complex ^, It has also been theorized that quark matter (known as strange quark matter) 
may be even more stable than atomic nuclei. In the latter event neutron stars should be entirely 
made of strange quark matter, possibly enveloped in a very thin nuclear crust. Such objects are 
called strange stars |||, ^ |lO|]. Strangeness, therefore, carried by hyperons, mesons, H-dibaryons, 
and strange quark matter, plays a key role for compact star physics and phenomenology, as will be 
discussed in this paper (see Ref. [Q] for a detailed recent review on this topic). 

2. Composition of high-density neutron star matter and EoS 

The properties of neutron stars are determined by the equation of state (EoS) of neutron star 
matter. The EoS of neutron star matter below neutron drip, which occurs at densities around 
4 X lO'^g/cm^, and at densities above neutron drip but below the saturation density of nuclear 
matter is relatively well known. This is to a less extent the case for the EoS in the vicinity of the 
saturation density of normal nuclear matter, = 0. 16 fm^^ (energy density of £0 = 140 MeV/fm'^). 
Finally, the physical properties of matter at still higher densities are extremely uncertain so that the 
associated EoS is only very poorly known ^, ||]. This is graphically indicated by the hatched 
areas in Fig. |lj, which is based on three Walecka-type models for the EoS. The HV model is a 
relativistic non-linear mean-field equation of state computed for nucleons (N= n, p) and hyperons 
(H= £, A, S) in chemical equilibrium [ll]]. The other two models, GfJo^" % and CFL [||], 
account additionally for the presence of up, down, and strange quarks in neutron star matter. The 
quarks are treated as normal unpaired quarks (Q) in the model, and as color-flavor locked 

superconducting quarks (sQ) with a superfluid gap of A = 100 MeV in the CFL model. The model 
labeled DD-RBHF is a density dependent relativistic Brueckner-Hartree-Fock EoS which accounts 
for nucleons and hyperons [|T3[ 0]. It is obvious from Fig. |I] that, depending on stellar composition, 
neutron star properties such as masses, radii, moments of inertia, redshifts, or limiting rotational 
periods may vary significantly with strangeness. 
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Figure 1: Models for the EoS of neutron star matter. 



2.1 Hyperons 

At the densities in the interior of neutron stars, the neutron chemical potential, /i", is likely 
to exceeds the masses, modified by interactions, of £, A and possibly E hyperons [11]. Hence, in 
addition to nucleons, neutron star matter is expected to have significant populations of hyperons 
and possibly even A's [15]. If so, pure neutron matter would constitute an excited state relative to 
hyperonic (many-baryon) matter which, therefore, would quickly transform via weak reactions like 



p + e +Ve 



(2.1) 



to the lower energy state. The chemical potentials associated with reaction ( p.lD in equilibrium 
obey the relation 

M" = A^^ + At'^ (2.2) 

where pL^' = since the mean free path of (anti) neutrinos is much smaller than the radius of 
neutron stars. Hence (anti) neutrinos do not accumulate inside neutron stars. This is different for 



hot proto-neutron stars []16|]. Equation ( |2.2[ ) is a special case of the general relation 

Ai^ =BV-^^Ai' , (2-3) 

which holds in any system characterized by two conserved charges. These are in the case of neutron 
star matter electric charge, q'"-, and baryon number charge, . Application of Eq. ( |23| ) to the A 
hyperon {B^ = \, = 0), for instance, leads to 

At^ = At" • (2.4) 
Ignoring particle interactions, the chemical potential of a relativistic particle of type % is given by 

Ai^ = (0{kF^) = ^ml + k\, (2.5) 

where (o{kf^ ) is the single-particle energy of the particle and kp^^ its Fermi momentum. Substituting 
( ^ ) into (^ leads to 

h 3 
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where = 1116 MeV and m„ = 939 MeV was used. That is, if interactions among the particles 
are ignored, neutrons are replaced with A's in neutron star matter at densities as low as two times the 
density of nuclear matter. This result is only slightly altered by the inclusion of particle interactions 
[jlTl]. Aside from chemical equilibrium, the condition of electric charge neutrality of neutron star 
matter, 

£ k\ + 3 0(/l^ - mM) = , (2.7) 

;t:=pX*,S",A++,...;e",^" 

where M stands for or K mesons, plays a key role for the particle composition of neutron 



star matter too. The last term in (2.7) accounts for the possible existence of either a ;r or a 



K meson condensate in neutron star matter, which will be discussed in more detail in Sect. 2.2 



below. Before, however, we illustrate the importance of Eqs. ( |2.2[ ) and ( |2.7[ ) for the proton-neutron 
fraction of neutron star matter. The beta decay and electron capture processes among nucleons, n —>■ 
p + e^ + Ve and p + e^ ^n + Vg respectively, also known as nucleon direct Urea processes, are only 
possible in neutron star matter if the proton fraction exceeds a certain critical value [ p^ . Otherwise 
energy and momentum can not be conserved simultaneously for these reactions so that they are 
forbidden. For a neutron star made up of only nucleons and electrons, it is rather straightforward to 
show that the critical proton fraction is around 11%. This follows from \ip^^ = k/r, + k/r, combined 
with the condition of electric charge neutrality of neutron star matter. The triangle inequality then 
requires for the magnitudes of the particle Fermi momenta ^f,, < kp^ +kf^, and charge neutrality 
dictates that kf^^ = kf^. Substituting kf^^ = kp^ into the triangle inequality leads to kp,, < Ikp^ so 
that for the particle number densities of neutrons and protons ?i„ < 8«p. Expressed as a fraction 
of the system's total baryon number density, n = rip +nn, one thus arrives at rip/n > 1/9 ~ 0.11, 
which is the figure quoted just above. Medium effects and interactions among the particles modify 
this value only slightly but the presence of muons raises it to about 0.15. Hyperons, which may 
exist in neutron star matter rather abundantly, produce neutrinos via direct Urea processes like 
^ A + e^ + Ve and A + e^ —>-'L^ + Ve [|l^] . The direct Urea processes are of key importance for 
neutron star cooling, which will be discussed briefly in Sect. ^ In most cases, the nucleon direct 



Urea process is more efficient than the ones involving hyperons [|19|, [20| ]. 

2.2 Meson condensates 

The pion or kaon meson fields may develop condensates in dense neutron star matter. These 
condensates would have two important effects on neutron stars. Firstly, they would soften the EoS 
above the critical density for onset of condensation, which reduces the maximum neutron star mass. 
Secondly, since the <n^> or <K > condensates can absorb as little or as much momentum as 
required by the scattering processes n+ < >^ n + e^ + Vg or n+ <K >^ n + + Vg, the 



associated neutrino emissivities are very high which leads to fast neutron star cooling [g0|, |21|] 
(see Sect. Since the K condensate process involves a change in strangeness, it is roughly 
by a factor sin^dc — 1/20 {dc denotes the Cabibbo angle) less efficient than the condensate 
process. However, medium effects can reduce the impact of the condensate on stellar cooling 
by about one order of magnitude, making it comparable to the efficiency of the K condensate. 
Estimates predict the onset of condensation at densities around ~ 2no, with no = 0.16 fm^^ 
the empirical nuclear matter density. However, this density is very sensitive to the strength of the 
effective nucleon particle-hole repulsion in the isospin T = I, spin S = I channel, which tends to 
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suppress tt condensation and may push to much higher values. Similarly, depending on the 
nuclear model, the threshold density for the onset of K condensation, , is probably at least as 



high as 4«o \9-A |23| ]. 

K condensation can only occur in neutron star matter if the electron chemical potential equals 
the effective in-medium meson mass, according to the schematic reaction K + v^, with the 
neutrinos leaving the star (see Fig. 0). This reaction would be followed by n ^ p + K^ . By this 
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Figure 2: The effective kaon mass, ot^, in nuclear matter and neutron star matter [Q]. Data taken from [ p4[ 
and [E2|, respectively. 



conversion the nucleons in the cores of newly formed neutron stars can become half neutrons and 
half protons [25|. The relatively isospin symmetric composition achieved in this way resembles 
the one of isospin symmetric atomic nuclei, which are made up of equal numbers of neutrons and 
protons. Neutron star matter is therefore referred to in this picture as nucleon matter, and neutron 
stars constructed for such an EoS are referred to as nucleon stars [25, 26, 27, 28]. 



2.3 H-dibaryons 

A novel particle that could make its appearance in the center of a neutron star is the so-called 
H-dibaryon, a doubly strange six-quark composite with spin and isospin zero, and baryon number 
two [p9|]. Since its first prediction in the 1970s, the H-dibaryon has been the subject of many 
theoretical and experimental studies as a possible candidate for a strongly bound exotic state. In 
neutron stars, which may contain a significant fraction of A hyperons, the A's could combine 
to form H-dibaryons, which could give way to the formation of H-dibaryon matter at densities 
somewhere above ~ 3 £o [|3^, 0, depending on the in-medium properties of the H-dibaryon. 
For an attractive optical potential, Uy{, of the H-dibaryon at normal nuclear density the equation 
of state is softened considerably, as shown in Fig. ^. H-dibaryon matter could thus exist in the 
cores of moderately dense neutron stars. H-dibaryons with a vacuum mass of about 2.2 GeV and 
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Energy Density [MeV fm ^] 

Figure 3: EoS of neutron star matter accounting for a H-dibaryon condensate |j3^ . Un is the optical potential 
of the H-dibaryon at normal nuclear density. 



a moderately attractive potential in the medium of about Un = —30 MeV, for instance, could go 
into a boson condensate in the cores of neutron stars if the limiting star mass is about that of the 
Hulse-Taylor pulsar PSR 1913+16, M = I.AAAMq [^]. Conversely, if the medium potential were 
moderately repulsive, around ?7h = +30 MeV, the formation of H-dibaryons may only take place in 
heavier neutron stars of masses greater than about \.6Mq. If formed, however, H-matter may not 
remain dormant in neutron stars but, because of its instability against compression, could trigger 
the conversion of neutron stars into hypothetical strange stars [31, 33, 



2.4 Quark deconflnement 

One item that came recently into particular focus concerns the possible existence of quark 
matter in the cores of neutron stars ||^ |5[ The phase transition from confined hadronic matter to 
deconfined quark matter is characterized by the conservation of baryon charge and electric charge. 
The Gibbs condition for phase equilibrium then is that the two associated chemical potentials, jU" 
and /l^, and the pressure in the two phases be equal |[l], I 



w^A^^m,^)=PQ(At^At^^). as) 

P]\ denoted the pressure of hadronic matter computed for a given hadronic matter Lagrangian 
^m{{x}), where {%} denotes the field variables and Fermi momenta that characterize a solution 
to the field equations of confined hadronic matter, 

{iy^ dp. - ) (^) = £ ^^Mx^ix) Yx (x) , (2.9) 
M=a,a),7i,... 

{d^d^+ml)a{x) = £ rax¥xi^)¥xix)^ (2-10) 

X=p,nX,- 

plus additional equations for the other meson fields (M = CO, 7r,p, ...). The pressure of quark matter, 
Pq, is obtainable from the bag model. The quark chemical potentials /x", jj,^, are related to the 
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r (km) r (km) r (km) 

Figure 4: Dependence of stellar compositions on neutron star spin frequency, V, for the HV, GgJ^", and CFL 
(from left to right) EoSs. The non-rotating stellar mass in each case is about \AMq. Vk denotes the Kepler 
(mass-shedding) frequency of each sequence, discussed in Sect. 0. 



baryon and charge chemical potentials as 



12 11 



(2.11) 



Equation ( |2.8[ ) is to be supplemented with the two global relations for conservation of baryon charge 
and electric charge within an unknown volume V containing A baryons. The first one is given by 



n = - = {\-^)ny,{pL\pL' ,T) + ^n^{pL\pL' ,T) , 



(2.12) 



where r] = Vq/V denotes the volume proportion of quark matter, Vq, in the unknown volume V , and 
?iH and HQ are the baryon number densities of hadronic matter and quark matter. Global neutrality 
of electric charge within the volume V can be written as 







Q 
V 



(2.13) 



with <7, the electric charge densities of hadrons, quarks, and leptons. For a given temperature, T, 
Eqs. ( ^^ to ( |2.13 ) serve to determine the two independent chemical potentials and the volume 
V for a specified volume fraction Tj of the quark phase in equilibrium with the hadronic phase. 
After completion Vq is obtained as Vq = rjV. Because of Eqs. (2^) through ( 2.13 ) the chemical 
potentials depend on the proportion rj of the phases in equilibrium, and hence so also all properties 
that depend on them, i.e. the energy densities, baryon and charge densities of each phase, and 
the common pressure. For the mixed phase, the volume proportion of quark matter varies from 
< 77 < 1 and the energy density is the linear combination of the two phases [[T], ^5[ ], 



£ = (l-r7)£H(M^M^w,n+T]£Q(At^A^^n• 



(2.14) 



Hypothetical neutron star compositions computed along the lines described above are shown in 
Fig. Possible astrophysical signals originating from quark deconfinement will be discussed in 
Sect.|§|[^,0]. 
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2.5 Color superconductivity of quark matter 

There has been much recent progress in our understanding of quark matter, culminating in the 
discovery that if quark matter exists it ought to be in a color superconducting state |6|, ^ 
This is made possible by the strong interaction among the quarks which is very attractive in some 
channels. Pairs of quarks are thus expected to form Cooper pairs very readily. Since pairs of quarks 
cannot be color neutral, the resulting condensate will break the local color symmetry and form what 
is called a color superconductor. The phase diagram of such matter is expected to be very complex 
[^, 1^. This is caused by the fact that quarks come in three different colors, different flavors, and 
different masses. Moreover, bulk matter is neutral with respect to both electric and color charge, 
and is in chemical equilibrium under the weak interaction processes that turn one quark flavor into 
another. To illustrate the condensation pattern briefly, we note the following pairing ansatz for the 
quark condensate [Q, 

{Wfcy^wl) ~ ^,£-^'efJ,,+^2e''^'\J,^ + ^,e-^hfJ,, , (2.15) 

where Yf is a quark of color a = {r,g,b) and flavor fa = {u,d,s), and Sjjk denotes the Levi-Civita 
symbol. The latter is zero for i = j,j = k, or k = i; +1 for {i,j,k) an even permutation of (1,2,3); 
and — 1 for (/, j, k) an odd permutation of (1,2,3). The condensate is a Lorentz scalar, antisymmetric 
in Dirac indices, antisymmetric in color, and thus antisymmetric in flavor. The gap parameters 
Ai, A2 and A3 describe d-s, u-s and u-d quark Cooper pairs, respectively. The following pairing 
schemes have emerged. At asymptotic densities (m^ ^ or jU ^ 00) the ground state of QCD with 
a vanishing strange quark mass is the color-flavor locked (CFL) phase (color-flavor locked quark 
pairing), in which all three quark flavors participate symmetrically. The gaps associated with this 
phase are 

A3~A2 = Ai=A, (2.16) 
and the quark condensates of the CFL phase are approximately of the form 

{Wlcyswl)-^£''^''efj„x, (2.17) 
with color and flavor indices all running from 1 to 3. Since e^l^^Ef^fi^x = ^/ ~ ^fb^f ^^^^ 



that the condensate (2.17) involves Kronecker delta functions that link color and flavor indices. 
Hence the notion color-flavor locking. The CFL phase has been shown to be electrically neutral 
without any need for electrons for a significant range of chemical potentials and strange quark 



masses [|41[]. If the strange quark mass is heavy enough to be ignored, then up and down quarks 
may pair in the two-flavor superconducting (2SC) phase. Other possible condensation patterns 
are CFL-i«:° CFL-K+ and CFL-tt" " gCFL (gapless CFL phase) ISC (single- 



flavor- pairing) ||4|,|4|, ||], CSL (color -spin locked phase) pq], and the LOFF (crystalline pairing) 



[g7|, |8|, ||] phase, depending on m.v, n, and electric charge density. Calculations performed for 
massless up and down quarks and a very heavy strange quark mass (m^ —>■ 00) agree that the quarks 
prefer to pair in the two-flavor superconducting (2SC) phase where 

A3>0, and A2=Ai=0. (2.18) 

In this case the pairing ansatz ( 2.151 ) reduces to 

{YfCYswi) °^ ASabE"'^' . (2.19) 
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Here the resulting condensate picks a color direction (3 or blue in the example ( |2.19 ) above), and 



creates a gap A at the Fermi surfaces of quarks with the other two out of three colors (red and 
green). The gapless CFL phase (gCFL) may prevail over the CFL and 2SC phases at intermediate 
values of m^//x with gaps given obeying the relation 

A3>A2>Ai>0. (2.20) 

For chemical potentials that are of astrophysical interest, /i < 1000 MeV, the gap is between 50 and 
100 MeV. The order of magnitude of this result agrees with calculations based on phenomenolog- 
ical effective interactions [ p9| , ^ ] as well as with perturbative calculations for /i > 10 GeV [pT|]. 
We also note that superconductivity modifies the equation of state at the order of (A/jj.)^ [ |T2| , ^], 
which is even for such large gaps only a few percent of the bulk energy. Such small effects may 
be safely neglected in present determinations of models for the equation of state of quark-hybrid 
stars. There has been much recent work on how color superconductivity in neutron stars could 
affect their properties (see Refs. ||5|,|6, 47,|53, 54, 55] and references therein). These studies reveal 



that possible signatures include the cooling by neutrino emission, the pattern of the arrival times 
of supernova neutrinos, the evolution of neutron star magnetic fields, rotational stellar instabilities, 
and glitches in rotation frequencies. 

Aside from neutron star properties, an additional test of color superconductivity may be pro- 
vided by upcoming cosmic ray space experiments such as AMS and ECCO [^].^ As shown 



in Ref. [58|, finite lumps of color-flavor locked strange quark matter (see Sect. 2.6), which should 
be present in cosmic rays if strange matter is the ground state of the strong interaction, turn out 
to be significantly more stable than strangelets without color-flavor locking for wide ranges of 
parameters. In addition, strangelets made of CFL strange matter obey a charge-mass relation of 
Z/A oc A^^/^, which differs significantly from the charge-mass relation of sti^angelets made of or- 
dinary strange quark matter. In the latter case, Z/A would be constant for small baryon numbers A 
and Z/A oc A^^/s for large A [|l^, ||, This difference may allow an experimental test of CFL 



locking in strange quark matter [58] 



2.6 Absolute stability of strange quark matter 

So far we have assumed that quark matter forms a state of matter higher in energy than atomic 
nuclei. This most plausible assumption, however, may not be correct [ |60| , |6T] , |62| ], since for a 
collection of more than a few hundred u, d, s quarks, the energy per baryon, E/A, of quark matter 
can be just as well below the energy per baryon of the most stable atomic nuclei, nickel and iron. 
This is known as the strange quark matter hypothesis. The energy per baryon in ^^Fe, for instance, 
is given by M(^^Fe)c^/56 = 930.4 MeV, with M(^^Fe) the mass of the ^^Fe atom. A simple 
estimate shows that for strange quark matter described by the MIT bag model E/A = ABn^/n^, 
so that bag constants of B = 57 MeV/fm^ (i.e. B^/"* = 145 MeV) and S = 85 MeV/fm^ (fii/^ = 
160 MeV) would place the energy per baryon of strange quark matter at E/A = 829 MeV and 
915 MeV, respectively, which correspond to strange quark matter which is absolutely bound with 
respect to nuclear matter [|l0|]. If this were indeed the case, neutron star matter would be metastable 
with respect to strange quark matter, and all neutron stars should in fact be strange quark stars 



See J. Madsen's contribution on strange matter in cosmic rays published elsewhere in this volume. 
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e„„„,=0.24 MeV/fm" i 




e (MeV/fm^) 



Figure 5: EoS of strange quark matter surrounded by nuclear matter. The maximal possible nuclear matter 
density is determined by neutron drip which occurs at Ecmst — 

0.24 MeV/fm^ (4.3 x lO" g/crci). Any 
nuclear density that is smaller than that is possible. As an example, we show here the EoS for a chosen 
density of gcmst = lO""* MeV/fm^ (lO'* g/cm^). 



[10, 63, 54|. As outlined just above, strange quark matter is expected to be a color superconductor 
which, at extremely high densities, should be in the CFL phase. This phase is rigorously electrically 
neutral with no electrons required [41]. For sufficiently large strange quark masses, however, the 
low density regime of strange quark matter is rather expected to form a 2SC phase (or possibly other 
phases) in which electrons are present [^, ^. The presence of electrons causes the formation of an 
electric dipole layer on the surface of strange matter, which enables strange quark matter stars to 
carry crusts made of ordinary nuclear matter [|8[ ^, |^]. The maximal possible density at the base of 
the crust (inner crust density) is determined by neutron drip, which occurs at about 4 x 10'' g/cm^. 
This somewhat complicated situation of the structure of strange matter enveloped in a (chemically 
equilibrated) nuclear crust can be represented by a proper choice for the EoS which consists of 
two parts []66|]. At densities below neutron drip it can be represented by the EoS of Baym-Pethick- 
Sutherland, while the high-density part, consisting of strange quark matter, can be described by 
the bag model EoS (see Fig. |5|). The EoS is characterized by a significant discontinuity in density 
between strange quark matter and nuclear crust matter across the electric dipole gap where the 
pressure of the nuclear crust at its base equals the pressure of strange matter at its surface ^ |6^ ]. 

One crucial astrophysical test of the strange quark matter hypothesis is whether strange quark 
stars can give rise to the observed phenomena of pulsar glitches. In the crust quake model an oblate 
solid nuclear crust in its present shape slowly comes out of equilibrium with the forces acting on it 
as the rotational period changes, and fractures when the built up stress exceeds the sheer strength 
of the crust material. The period and rate of change of period slowly heal to the trend preceding 
the glitch as the coupling between crust and core re-establish their co-rotation. The existence of 
glitches may have a decisive impact on the question of whether the strange quark matter hypothesis 
holds or not. From the calculations in [|6^] it is known that the ratio of the crustal moment of inertia 
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to the star's total moment of inertia, /cmst/Aotah varies between 10^^ and 10^^ at the maximum 
mass. If the angular momentum of the pulsar is conserved in the quake then the relative frequency 
change and moment of inertia change are equal and one arrives at 



Ail |A/| |A/| /(-riisf . <; 5^ 

— = ^ > ^ ^ (10-5-10-3)/, with 0</<l. (2.21) 

Here Iq denotes the moment of inertia of that part of the star whose frequency is changed in the 
quake. It might be that of the crust only, some fraction, or all of the star. The factor / in Eq. 
( |2.21 ) represents the fraction of the crustal moment of inertia that is altered in the quake, i.e., / = 



|A/|//cmst- Since the observed glitches have relative frequency changes = (10 ^ — 10 ^), 

a change in the crustal moment of inertia of / 5 0.1 would cause a giant glitch even in the least 



favorable case |66]. Moreover, one finds that the observed range of the fractional change in the 



spin-down rate, Cl, is consistent with the crust having the small moment of inertia calculated and 



the quake involving only a small fraction / of that, just as in Eq. (2.21). To this aim we write JSq] 



\Q. _ \Q./Q. \M\ _ \Q./Q. ^ /crust 



Cl ^Q./Q. Iq AQ./Q. 



^^rusi ^ (10-1 to 10)/, (2.22) 



where use of Eq. ( 2.21 ) has been made. Equation ( 2.22 ) yields a small / value in the range / < 
(10-"* to 10-^), in agreement with / 5 10-' established just above. Here measured values of the 
ratio {AD./Q.) / {AQ./Q.) ~ 10-^ to 10—* for the Crab and Vela pulsars, respectively, have been used. 
An improved discussion of the surface gap below strange star crusts has been performed very 



recently in Ref. []67p. In addition to the electrostatic forces described above, this study includes 
gravity too. The properties of the gap are investigated for a wide range of parameters assuming both 
color-flavor locked and noncolor-flavor locked strange star cores. It is found that the maximally 
allowed inner crust density is generally lower than that of neutron drip. This does not alter the 
overall form of the EoS shown in Fig. ^, however. Another interesting finding concerning the 
surface properties of strange stars was recently published in Ref. [|68|]. In this paper it is found that, 
depending on the surface tension of nuggets of strange matter, a heterogeneous crust comprised 
of nuggets of strange quark matter embedded in an uniform electron background may exist in the 
surface region of strange stars. This heterogeneous strange star surface would have a negligible 
electric field which would make the existence of an ordinary nuclear crust, which requires a very 
strong electric field, impossible. 



3. Models of compact stars 

Neutron stars are objects of highly compressed matter so that the geometry of space-time 
is changed considerably from flat space. Thus models of such stars are to be constructed in the 
framework of Einstein's general theory of relativity combined with theories of superdense matter. 
The effects of curved space-time are included by coupling the energy-momentum density tensor 
for matter fields to Einstein's field equations. The generally covariant Lagrangian density is 

^ = ^E + ^G, (3.1) 
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where the dynamics of particles is introduced through the matter Lagrangian added to the 
gravitational Lagrangian ^q. The latter is given by 

^G=g"'R = g"h^''Rl.y, (3.2) 

where g^^ and R^y denote the metric tensor and the Ricci tensor, respectively. The latter is given 
by 

where commas followed by a Greek letter denote derivatives with respect to space-time coordinates, 



e.g. ,v = d /dx^ etc. The Christoffel symbols T in (3.3) are defined as 



r;iv = (gM,v +8vX^n-8^v,x) . (3.4) 

The connection between both branches of physics is provided by Einstein's field equations 

= Rtiv _ Igi^yR = SKT^\e,P{e)) , (3.5) 

(/X, V = 0, 1,2,3) which couples the Einstein curvature tensor, G'^^, to the energy-momentum den- 



sity tensor, T^^, of the stellar matter. The quantities g^"^ and R in ( |3.5| ) denote the metric tensor 
and the Ricci scalar (scalar curvature) The tensor T^^ contains the equation of state, P{£), 
of the stellar matter discussed in Sect. |^ In general, Einstein's field equations and the many-body 
equations were to be solved simultaneously since the baryons and quarks move in curved space- 
time whose geometry, determined by Einstein's field equations, is coupled to the total mass energy 
density of the matter. In the case of neutron stars, as for all astrophysical situations for which the 
long-range gravitational forces can be cleanly separated from the short-range forces, the deviation 
from flat space-time over the length scale of the strong interaction, ~ 1 fm, is however practically 
zero up to the highest densities reached in the cores of such stars (some 10^^ g/cm^). This is not 
to be confused with the global length scale of neutron stars, ~ 10 km, for which M/R ~ 0.3, de- 
pending on the star's mass. That is to say, gravity curves space-time only on a macroscopic length 
scale but leaves it flat to a very good approximation on a microscopic length scale. To achieve an 
appreciable curvature on a microscopic scale set by the strong interaction, mass densities greater 



than ~ lO^'' g/cm would be necessary [69]! This circumstance divides the construction of models 
of compact stars into two distinct problems. Firstly, the effects of the short-range nuclear forces on 
the properties of matter are described in a comoving proper reference frame (local inertial frame), 
where space-time is flat, by the parameters and laws of special relativistic many-body physics. 
Secondly, the coupling between the long-range gravitational field and the matter is then taken into 
account by solving Einstein's field equations for the gravitational field described by the general rel- 
ativistic curvature of space-time, which determines the global structure of the stellar configuration. 

For many studies of neutron star properties it is sufficient to treat neutron star matter as a 
perfect fluid. The energy-momentum tensor of such a fluid is given by 

T^^ = ^^{£ + p)+g^^P. (3.6) 
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Figure 6: Mass-radius relations of non-rotating (left panel) and rotating (right panel: fi = Q.k) neutron stars 
computed for different EoSs. £Ik denotes the general relativistic Kepler (mass shedding) frequency. 



For non-rotating spherically symmetric stars the metric has the rather simple form 

ds^ = -e^'^t'') dt'^ + dr^ + dd^ + sin^Q , (3.7) 



where <I>(r) and A(r) are radially varying metric functions. From ( |3.7[ ) one reads off the following 
covariant components of the metric tensor, 

ga = -e^""^'^ , grr = e'""^'^ , gee = , g^ = '-'sin'^ , (3.8) 

so that the only non- vanishing Christoffel symbols are 

n = ^'{r) , rl = ^'{r) , F^, = A'(r) , r% = r% = r'^ , F^,, = - r .-2A(^) , 

= ^ , = - r sin^e ^^-2A(,) , = - sin COS , (3.9) 

where primes denote differentiation with respect to the radial coordinate. Substituting (|3.6|) and 



(3.S) into Einstein's field equations leads to the general relativistic equations of hydrostatic equi- 



(3.10) 



librium discussed first by Tolman ||70p and Oppenheimer-Volkoff [|71|], 

dP{r) _ e{r)m{r) y + Wi) V '"W 



dr f- \_ 

r 

Note that we use geometrized units, where the gravitational constant and velocity of light are 
G = c = 1 so that Mq = 1 .475 km. The boundary condition to ( |3.10[ ) is P{r = 0) = P(£c), where e,. 
denotes the energy density at the star's center, which constitutes a parameter. Equation ( [3.10 ) is to 
be integrated out to a radial distance where P(r) = which determines the star's radius, R. The mass 
contained in a sphere of radius r (< R), denoted by m(r), follows from m(r) = 4?! Jgd/ r'^ £(/) . 
The star's total gravitational mass is thus given by M = m{R). Figure ^ shows the mass-radius 
relationships of both non-rotating as well as rotating sequences of neutron stars for the sample 
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EoSs discussed in Sect. |2[ The non-rotating sequences are solutions of the Tolman-Oppenheimer- 
Volkoff equation shown in ( |3.10| ). The construction of rotating sequences will be discussed shortly 
below. Figure ^ shows the gravitational mass of non-rotating as well as rotating neutron stars as 
a function of central star density. Stars to the right of the respective mass peaks in each panel are 
unstable against radial oscillations and thus cannot exist stably in nature. Also shown in these plots 
are the evolutionary (constant stellar baryon number, A) paths that isolated rotating neutron stars 
would follow during their stellar spin-down evolution caused by the emission of magnetic dipole 
radiation and a wind of pairs. Figure ^ reveals that CFL stars may spend considerably more 

time in the spin-down phase than their competitors of the same mass. Another point that we want to 




100 1000 100 1000 100 1000 

e (MeV/fm^) 



Figure 7: Mass-central energy relations for the sample equations of state introduced in Sect. ^. 



make is that all equations of state are able to support neutron stars of canonical mass, M ~ IAMq. 
Neutron stars more massive than about IMq,, on the other hand, are only supported by equations of 
state that exhibit a very stiff behavior at asymptotic densities, disfavoring the presence of hyperons, 
meson condensates, or quarks. Knowledge of the maximum possible mass of neutron stars is of 
great importance for two reasons. Firstly, because the largest known neutron star mass imposes a 
lower bound on the maximum mass of a theoretical model. Very massive neutron star candidates 
are J0751-F 1807 (I.I^^^Mq @]), VelaX-1 (l.SSiO.lSM© if the inclination angle of the system 
is / = 90°; an inclination angle of / = 70° increases the star's mass to 2.27 ztO-HMg, [[73|]), and 
Cyg X-2 (1.78 ±0.23 Mq [74]. Titarchuck and Shaposhnikov obtain for Cyg X-2 a lower mass of 
1.44ibO.O6M0 [75 1). The second reason is that the maximum mass of neutron stars is essential in 
order to identify solar-mass black hole candidates [76, 77 1. 

The structure equations of rotating compact stars are considerably more complicated that those 
of non-rotating compact stars [^. These complications have their cause in the rotational deforma- 
tion, that is, a flattening at the pole accompanied with a radial blowup in the equatorial direction, 
which leads to a dependence of the star's metric on the polar coordinate, 6, in addition to the 
mere dependence on the radial coordinate, r. Secondly, rotation stabilizes a star against gravi- 
tational collapse. A rotating star can therefore carry more mass than a non-rotating star. Being 
more massive, however, means that the geometry of space-time is changed too. This makes the 
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Figure 8: Dragging of local inertial frames (Lense-Thirring effect) caused by 1 .4M0 neutron stars rotating 



at 2 ms. The frequency (O is defined in Eq. (3.13) 



metric functions associated with a rotating star depend on tiie star's rotational frequency. Finally, 
the general relativistic effect of the dragging of local inertial frames implies the occurrence of an 
additional non-diagonal term, , in the metric tensor g^^ . This term imposes a self-consistency 
condition on the stellar structure equations, since the extent to which the local inertial frames are 
dragged along by the star is determined by the initially unknown stellar properties like mass and 
limiting rotational frequency. The covariant components of the metric tensor of a rotating compact 



star are thus given by [E, 78] 



gtt 



„2A 



,2 V 



which leads for the line element to 



ds^ = g^ydx^dx" = -e^''dt'- + e^'^ {d^ -(odtf + e^^ dd^ + e^^ dr^ . 



(3.11) 



(3.12) 



Here each metric function, i.e. v, i/A, /i and A, as well as the angular velocities of the local inertial 
frames, (O, depend on the radial coordinate r and polar angle Q and implicitly on the star's angular 
velocity Q.. Of particular interest is the relative angular frame dragging frequency, d), defined as 



(y(r,0,a) =a-w(r,0,a), 



(3.13) 



which is the angular velocity of the star, Q., relative to the angular velocity of a local inertial frame, 
ft). It is this frequency that is of relevance when discussing the rotational flow of the fluid inside the 
star, since the magnitude of the centrifugal force acting on a fluid element is governed-in general 
relativity as well as in Newtonian gravitational theory-by the rate of rotation of the fluid element 
relative to a local inertial frame [[7^]. In contrast to Newtonian theory, however, the inertial frames 
inside (and outside) a general relativistic fluid are not at rest with respect to the distant stars, as 
pointed out just above. Rather, the local inertial frames are dragged along by the rotating fluid. 
Depending on the internal stellar constitution, this effect can be quite strong, as shown in Fig. |8]for 
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Table 1: Properties of neutron stars composed of nucleons and hyperons (HV), nucleons, hyperons, and 
normal quarks (Ofgo")' ™d nucleons, hyperons, and color-superconducting quarks (CFL). 
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§8,14 (cm/s2) 
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3.0146 


0.7278 
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1.4493 


T/W 











0.0894 


0.0941 


0.0787 


BE (Me) 


0.0937 


0.1470 


0.1534 


0.0524 


0.1097 


0.1203 













0.336 


0.353 


0.424 



rotating 2 ms neutron stars. For a very compact neutron star containing a color-superconducting 
CFL core, as in our example, one reads off from this figure that the local inertial frames at the star's 
center rotate at about half the star's rotational frequency, (o{r = 0) ~ Q./2. This value drops to 
about 15% for the local inertial frames located at the star's equator. The scenarios shown in Fig. 
|8| may be of great importance for binary millisecond neutron stars in their final accretion stages, 
where the accretion disk approaches the star very closely. 

Table [T] summarizes the impact of strangeness on several intriguing properties of non-rotating 
as well as rotating neutron stars. The latter spin at their respective Kepler frequencies. One sees 
that the central energy density, Ec, spans a very wide range, depending on particle composition. 
The surface redshift is of importance since it is connected to observed neutron star temperatures 
through the relation T°° /T^^f = 1/(1 +Z). CFL quark stars may have redshifts that are up to 50% 
higher than those of conventional stars. Finally, we also show in Table |l| the surface gravity of 
stars, gs.i4 [80], which again may be up to 50% higher for CFL stars. The other quantities listed are 
the rotational kinetic energy in units of the total energy of the star, T /W , the stellar binding energy, 
BE, and the rotational velocity of a particle at the star's equator 



4. Limiting rotational periods 

4.1 Mass shedding from the equator 

No simple stability criteria are known for rapidly rotating stellar configurations in general 
relativity. However, an absolute limit on rapid rotation is set by the onset of mass shedding from 
the equator of a rotating star. The corresponding rotational frequency is known as the Kepler 
frequency, Hk- In classical mechanics, the expression for the Kepler frequency, determined by the 
equality between the centrifugal force and gravity, is readily obtained as Hk = \/M/R^. In order 
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Figure 9: Moment of inertia of several sample neutron stars. 



to derive the general relativistic counterpart of this relation, one applies the extremal principle to 
the circular orbit of a point mass rotating at the star's equator. Since r = 6 = const for a point mass 



there, one has dr = d6 =0. The line element ( 3.12 ) then reduces to 

„2 



„2v 



(4.1) 



Substituting this expression into J = J^^- ds, where and S2 refer to points located at that particular 
orbit for which J becomes extremal, gives 



dt\ e^" -e^l'ia-cof. 



(4.2) 



Applying the extremal condition 57 = to Eq. ( |4.2[ ) and noticing that V = e^ ^ — co) then leads 
to the following relation. 



dco 



dv 



or or or 



(4.3) 



It constitutes a simple quadratic equation for the orbital velocity V of a particle at the star's equator. 
One thus obtains for the Kepler frequency (Kepler period, Pk) the final relation [0], 



(0' 



^K = (0 + — + e 



CO' 



V' ( „ 



1% 



(4.4) 



which is to be determined self-consistently at the star's equator (primes denote radial derivatives). 
For most neutron star matter equations of state, the Kepler period obtained for \AMq neutron stars 
scatters around 1 ms. One exception to this are strange quark matter stars. These are self-bound 
and, thus, tend to possess smaller radii than conventional neutron stars, which are bound by gravity 
only. Because of their smaller radii, strange stars can withstand mass shedding down to periods of 



around 0.5 ms [66, 81]. CFL stars reside between these limits. 
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Figure 10: Temperature dependence of the critical angular velocity €1^ of rotating neutron stars. The left 
panel shows the gravitational radiation driven /-mode instability suppressed by shear and bulk viscosity. 
Right panel; comparison of /-mode instabihty with r-mode instability. (Data from Refs. Js^, |83[].) 



As a last topic of this section, we briefly discuss the moment of inertia of a rotationally de- 



im =2nrder dr Z+^'+^+v' '-±^ ^ (4 5) 

^ ' Jo Jo ^'2v-2V-(«-n)2 a ■ 



formed star described by the metric in Eq. ( p.l2[ ). For such stars the moment of inertia is given 
by 

rTl rR(t 

Figure ^ shows that the crustal fraction of the moment of inertia of a neutron star may be around 
50% smaller if the star contains a very soft phase of matter like quark matter. This may be of 
relevance for pulsar glitch models and the modeling of the post-glitch behavior of pulsars. 

4.2 Gravitational radiation reaction driven instabilities 

Rotational instabilities in rotating stars, known as gravitational radiation driven instabilities, 
are probably setting a more stringent limit on rapid stellar rotation than mass shedding. These 
instabilities originate from counter-rotating surface vibrational modes which at sufficiently high 
rotational star frequencies are dragged forward. In this case gravitational radiation, which in- 
evitably accompanies the aspherical transport of matter, does not damp the instability modes but 
rather drives them. Viscosity plays the important role of damping these instabilities at a sufficiently 
reduced rotational frequency such that the viscous damping rate and power in gravity waves are 
comparable. The most critical instability modes that are driven unstable by gravitational radiation 
are /-modes and r-modes. Figure |T^ shows the stable neutron star frequencies if only /-modes 
were operative. One sees that hot as well as cold neutron stars can rotate at frequencies close to 
mass shedding, because of the large contributions of shear and bulk viscosity, respectively, for this 
temperature regime. The more recently discovered r-mode instability may change the picture com- 
pletely, as can be seen too from Fig. [T^. These modes are driven unstable by gravitational radiation 
over a considerably wider range of angular velocities than the /-modes (cf. dashed curve labeled 
(m = 2) r-mode instability). In stars with cores cooler than ~ 10^ K, on the other hand, the r-mode 
instability may be completely suppressed by viscous phenomena so that stable rotation would be 
limited by the /-mode instability again [Ml]. 
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Figure 11: Critical rotation frequencies versus 
stellar temperature for CFL strange stars [^5|]. 



Figure 12: Same as Fig. HJ but for 2SC quark 
stars [S. 



Figures ^ and 12 are the counterparts to Fig. 10 but calculated for strange stars made of 
CFL and 2SC quark matter, respectively [p4[ 85]. The r-mode instability seems to rule out that 



pulsars are CFL strange stars, if the characteristic time scale for viscous damping of r-modes are 



exponentially increased by factors of A/r as calculated in [84]. An energy gap as small as A = 

I MeV was assumed. For much larger gaps of A ~ 100 MeV, as expected for color superconducting 
quark matter in the CFL phase, the entire diagram would be r-mode unstable. The full curve in Fig. 

II is calculated for a strange quark mass of = 200 MeV, the dotted curve for = 100 MeV. The 



box marks the positions of most low mass X-ray binaries (LMXBs) []86|], and the crosses denote 
the most rapidly rotating millisecond pulsars known. All strange stars above the curves would 
spin down on a time scale of hours due to the r-mode instability, in complete contradiction to the 
observation of millisecond pulsars and LMXBs, which would rule out CFL quark matter in strange 
stars (see, however, Ref. [p7|]). Figure 12 shows the critical rotation frequencies of quark stars as 
a function of internal stellar temperature for 2SC quark stars. For such quark stars the situation is 
less conclusive. Rapid spin-down, driven by the r-mode gravitational radiation instability, would 
happen for stars above the curves. 



5. Astrophysical signals of quark deconfinement 

5.1 Isolated, rotating neutron stars 

Whether or not quark deconfinement occurs in neutron stars makes only very little difference 
to their static properties, such as the range of possible masses and radii, which renders the detection 
of quark matter in such objects extremely complicated. This may be strikingly different for isolated, 
rotating neutron stars which spin down, and thus become more compressed, because of the emission 
of magnetic dipole radiation and a wind of electron-positron pairs. For some rotating neutron stars 
the mass and initial rotational frequency may be just such that the central density rises from below 
to above the critical density for dissolution of baryons into their quark constituents. If accompanied 
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by a pronounced shrinkage of the neutron star, as is the case for the neutron star shown in the 
central panel in Fig. the star's moment of inertia could change dramatically. As shown in [^], 
the moment of inertia can decrease so anomalously that it could even introduce an era stellar spin- 
up that may last for ~ 10^ years. Since the dipole age of millisecond pulsars is about 10^ years, 
one may roughly estimate that about 10% of the solitary millisecond pulsars could be in the quark 
transition epoch and thus could be signaling the ongoing process of quark deconfinement. Changes 
in the moment of inertia reflect themselves in the braking index, n, of a rotating neutron star, as can 
be seen from [i, pi |89l] 



aa i + 3ra+i"a^ 



3rQ. + I"Q? 
21 + I'Q. 



(5.1) 



where dots (primes) denote derivatives with respect to time (D.). The last relation in ( |5.1| ) constitutes 
the non-relativistic limit of the braking index [|90|]. It is obvious that these expressions reduce to 
the canonical limit « = 3 if / is independent of frequency. Evidently, this is not the case for rapidly 



Table 2: Dominant neutrino emitting processes in neutron star cores if hyperons and quarks are absent [|21| 



Name 



Process 



Emissivity 
(erg cm^^ s^^) 



Modified Urea cycle 

(neutron branch) 
Modified Urea cycle 

(proton branch) 
Bremsstrahlung 

Cooper pair 
formations 
Direct Urea cycle 

condensate 
K condensate 



n + n ^ n + p + e +Vg 
n + p + e^ n + n + Ve 

p + n^ p_+ p + e^ + % 
p + p + e ^ p + n + Ve 

n + n—>n + n + v + v_ 
n+p^n+p+y+v 
p + p^ p + p + v + v 

n + n^ \nn] + V + V_ 
P + P^[pp] + V + V 

n p + e^ + Ve 
p + e ^n + Ve 



n+ <n >- 
n+ < K >- 



n + e + Vg 



2x10^^ RT^ 



~ ioi^/?r| 

^ 5xlO^\ RTl 
-5x10'^/? 77 

~ 10^'^ RT^ 

~ 10^^ RT^ 
~ 10^5 R 



Slow 
Slow 

Slow 

Medium 

Fast 

Fast 
Fast 



rotating neutron stars, and it also fails for stars that undergo pronounced compositional changes 
(phase transitions) which alter the moment of inertia significantly. Under favorable circumstances, 
these changes in I, originating from the transition of confined hadronic matter into quark matter, 
may cause the braking index to deviate dramatically from 3 in the vicinity of the star's frequency 
where the phase transition to quark matter occurs. The changes in I may even be so pronounced that 
n{D.) —5- ±00 at the transition point ||l|,|^, ^,^]. Such dramatic anomalies in n{Q.) are not known for 
conventional neutron stars (see left panel in Fig. because their radii and thus moments of inertia 
appear to vary smoothly with 11 [^, A counterexample to this, however, is discussed in [pT|]. 
The future astrophysical observation of strong anomalies in the braking behavior of isolated pulsars 
could thus be cautiously interpreted as a possible astrophysical signal for quark deconfinement in 
neutron stars. 
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Table 3: Dominant neutrino emitting processes in deconfined quark matter [Elh. 



Name 



Process 



Emissivity Efficiency 



(erg cm ^ s ^) 



Direct Urea cycle 

{ud branch) 
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{us branch) 
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dd 
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~ 10^''/? r| 

~ 1.5x102 /?7g 
~ 1.5x10^1 /?7^ 



Fast 

Fast 
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Slow 
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5.2 Accreting neutron stars 

Accreting x-ray neutron stars provide a very interesting contrast to the spin-down of isolated 
neutron stars. These x-ray neutron stars are being spun up by the accretion of matter from a lower- 
mass {M^OAMq), less-dense companion. If the critical deconfinement density falls within that of 
the canonical pulsars, quark matter could already exist in them but will be spun out of such stars 



as their frequency increases during accretion. This scenario has been modeled in [92], where it 
was found that quark matter remains relatively dormant in the core of a neutron star until the star 
has been spun up to frequencies at which the central density is about to drop below the threshold 
density at which quark matter exists. As known from the discussion above, this could manifest 
itself in a significant increase of the star's moment of inertia. The angular momentum added to a 
neutron star during this phase of evolution is therefore consumed by the star's expansion, inhibiting 
a further spin-up until the star's quark matter content has been completely converted into a mixed 
phase of hadrons and quarks. Such accreters, therefore, tend to spend a greater length of time in the 
critical frequencies than otherwise. For canonical accretion rates of 10^^° Mo/year the time span 
can be on the order of 10^ years. Hence, from this scenario, one would expect a greater number of 
accreting x-ray neutron stars that appear near the same frequency. Evidence that accreting neutron 
stars pile up at certain frequencies, which are well below the mass shedding limit, is provided by 
the spin distribution of accreting millisecond pulsars in 57 Tuc and neutron stars in low mass X-ray 
binaries observed with the Rossi X-ray Timing Explorer (RXTE). The proposed limiting mecha- 
nisms responsible for this behavior is generally attributed to gravity-wave emission caused by the 
r-mode instability, or by a small stellar mass quadrupole moment [93, 0, Supplemental to 
these explanations, quark reconfinement (or, more generally, strong first-order like phase transition) 
may be linked to this phenomenon as well [37, 92, 9^, p7| ]. 



6. Cooling of neutron stars 

The predominant cooling mechanism of hot (temperatures of several ~ lO'^ k) newly formed 
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Figure 13: Cooling behavior of a \AMq neutron star for competing assumptions about the properties of 
superdense matter. Three distinct cooling scenarios, referred to as standard, intermediate, and enhanced can 
be distinguished. The band-like structures reflect the uncertainties inherent in the stellar EoS Q|. 



neutron stars immediately after formation is neutrino emission, with an initial cooling time scale of 
seconds. Already a few minutes after birth, the internal neutron star temperature drops to ~ 10^ K. 
Photon emission overtakes neutrino emission when the internal temperature has fallen to ~ 10^ K, 
with a corresponding surface temperature roughly two orders of magnitude smaller. Neutrino cool- 
ing dominates for at least the first 10^ years, and typically for much longer in standard cooling 
(modified Urea) calculations. The dominant neutrino emitting processes in neutron star matter are 
summarized in Tables ^ and |[ Figure 13 shows the outcome of cooling calculations performed 
for a broad collection of equations of state [Q, ^ and competing assumptions about the dominant 
neutrino emitting processes. For recent overviews of neutron star cooling, see, for instance, Refs. 
ni |98ll. We also refer to D. Blaschke's contribution contained elsewhere in this volume. 



7. Summary 

It is often stressed that there has never been a more exciting time in the overlapping areas 
of nuclear physics, particle physics and relativistic astrophysics than today. This comes at a time 
where new orbiting observatories such as the Hubble Space Telescope (HST), Rossi X-ray Timing 
Explorer, Chandra X-ray satellite, and the X-ray Multi Mirror Mission (XMM) have extended 
our vision tremendously, allowing us to observe compact star phenomena with an unprecedented 
clarity and angular resolution that previously were only imagined. On the Earth, radio telescopes 
(Arecibo, Green Bank, Parkes, VLA) and instruments using adaptive optics and other revolutionary 
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techniques have exceeded previous expectations of what can be accompUshed from the ground. 
Finally, the gravitational wave detectors LIGO, LISA, and VIRGO are opening up a window for 
the detection of gravitational waves emitted from compact stellar objects such as neutron stars and 
black holes. This unprecedented situation is providing us with key information on compact stars, 
which are the only physical objects in which cold and dense baryonic matter is realized in nature. 
As discussed in this paper, a key role in compact star physics is played by strangeness. It alters 
the masses, radii, cooling behavior, and surface composition of neutron stars. Other important 
observables may be the spin evolution of isolated neutron stars and neutron stars in low-mass x-ray 
binaries. All told, these observables are key in exploring the phase diagram of dense nuclear matter 
at high baryon number density but low temperature, which is not accessible to relativistic heavy 
ion collision experiments. 
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